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Abstract 

In this paper we study point transformed electromagnetic invisibility cloaks in transformation 
media that are obtained by transformation from general anisotropic media. We assume that 
there are several point transformed electromagnetic cloaks located in different points in 
space. Our results apply in particular to the first order invisibility cloaks introduced by 
Pendry et al. and to the high order invisibility cloaks introduced by Hendi et al. and 
by Cai et al.. We identify the appropriate cloaking boundary conditions that the solutions 
of Maxwell equations have to satisfy at the outside, dK+, and at the inside, dK^, of the 
boundary of the cloaked object K in the case where the permittivity and the permeability 
are bounded below and above in K. Namely, that the tangential components of the electric 
and the magnetic fields have to vanish at dK^ - what is always true- and that the normal 
components of the curl of the electric and the magnetic fields have to vanish at dK_. These 
results are proven requiring that energy be conserved. In the case of one spherical cloak with 
a spherically stratified K and a radial current at dK we verify by an explicit calculation that 
our cloaking boundary conditions are satisfied and that cloaking of active devices holds, even 
if the current is at the boundary of the cloaked object. As we prove our results for media 
that are obtained by transformation from general anisotropic media, our results apply to the 
cloaking of objects with passive and active devices contained in general anisotropic media, 
in particular to objects with passive and active devices contained inside general crystals. 
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Our results suggest a method to enhance cloaking in the approximate transformation media 
that are used in practice. Namely, to coat the boundary of the cloaked object (the inner 
boundary of the cloak) by a material that imposes the boundary conditions above. As these 
boundary conditions have to be satisfied for exact transformation media, adding a lining that 
enforces them in the case of approximate transformation media will improve the performance 
of approximate cloaks. 

1 Introduction 

In this paper we study point transformed electromagnetic invisibility cloaks in transformation 
media that are obtained by transformation from general anisotropic media. We assume that 
there are several cloaks located in different points in space. Our results apply in particular 
to the first order invisibility cloaks introduced by [1] and to the high order invisibility cloaks 
introduced by [21 [3]. 

In [5] we gave a rigorous proof -based in energy conservation- of cloaking of pas- 
sive and active devices for our general class of invisibility cloaks. The cloaked object, 
completely decouples from the exterior. Actually, the cloaking outside is independent of 
what is inside K. The electromagnetic waves inside K can not leave K and vice versa, the 
electromagnetic waves outside can not go inside. Furthermore, we identified the appropriate 
cloaking boundary conditions when cloaking is formulated as a boundary value problem. We 
proved that the tangential components of the electric and the magnetic fields have to vanish 
at the outside of the boundary of the cloaked object, dK^. This boundary condition is self- 
adjoint in our case because the permittivity and the permeability are degenerate at dK^. 
We also proved that the boundary condition at the inside of the boundary of the cloaked 
object, dK_, can be any self-adjoint boundary condition for the Maxwell generator in K . 
This is true in the general case where the permittivity and the permeability are allowed to be 
degenerate at dK^ . In this general situation the particular boundary condition that nature 
will take on dK_ will depend on the behaviour of the permittivity and the permeability 
near dK_. We proved our results both in the time and in the frequency domains, and as 
we consider media obtained by transformation from general anisotropic media, our results 
apply, in particular, to objects inside general crystals. 

In this paper we address the problem of determining the cloaking boundary conditions at 
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dK_ in the case where the permittivity and the permeabihty inside K are bounded and have 
a positive lower bound. This corresponds to the situation where we have a standard object K 
-that could be anisotropic and inhomogeneous, but whose permittivity and permeability are 
neither singular nor degenerate- that is coated by a transformation medium that is degenerate 
at dKj^. We also allow for active devices in K. This is perhaps the more important case in 
the applications. We prove that in this case the cloaking boundary conditions at dK_ are 
that the normal components of the curl of the electric and the magnetic fields vanish. 

Since we have identified the cloaking boundary conditions at dK± we have now a complete 
formulation of cloaking as a boundary value problem in this important case. For the exact 
transformation media that we consider in this paper these boundary conditions are satisfied 
because they follow from energy conservation and there is no need to add any lining to impose 
them. In other words, they are the conditions taken by nature, as they are imposed by energy 
conservation. However, our results suggest a method to enhance cloaking in the approximate 
transformation media that are used in practice. Namely, to coat the boundary of the cloaked 
object (the inner boundary of the cloak) by a material that imposes the boundary conditions 
above. As these boundary conditions have to be satisfied for exact transformation media, 
adding a lining that enforces them in the case of approximate transformation media will 
improve the performance of approximate cloaks. 

It is, of course, a well known fact in electromagnetic theory -and in wave propagation 
in general- that in any interphase between two different media there has to be a boundary 
condition. This obviously applies to the interphase between the cloaked object and the 
coating metamaterial, i.e., dK. So, the real question is not if there has to be boundary 
conditions at dK±, but rather what are the appropriate boundary conditions. We address 
this question in this paper as well as in [HE]. The reason why this is a delicate problem 
that requires a careful analysis is that for point transformed electromagnetic cloaks the 
permittivity and the permeability are degenerate at dK and in consequence the standard 
rules that are used in the non-degenerate case do not apply. In fact, the solutions to Maxwell 
equations are, in general, discontinuous at dK. 

The interesting paper [6] considers point and line transformed electromagnetic cloaks 
under general coordinate transformations. Among other problems, they compute the fields 
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outside the cloaked object from the fields in the original electromagnetic space, using the 
transformation formulae between them, what avoids doing tedious calculations in the trans- 
formed space. This method was previously used in [H |3] for spherical and cylindrical cloaks. 
In this way it is proven in [B] that for general point transformed cloaks the tangential com- 
ponents of the electric and the magnetic fields vanish at the outside of the boundary of the 
cloaked object and also that for general line transformed cloaks the tangential components 
of the electric and magnetic fields that are orthogonal to the axis of the cloak vanish at the 
outside of the boundary of the cloaked object. This generalizes the results previously proved 
in [H [5] in the case of spherical and cylindrical cloaks, using the same method. 

The paper [7| considers cloaking in terms of the Cauchy data, in the context of the 
Dirichlet to Neumann operator. Among other problems, they study cloaking of passive 
and active devices for one spherical electromagnetic cloak. They postulate a class of weak 
solutions in distribution sense across dK (see Definition 4 of [7]). They study the case when 
the permittivity and permeability are bounded below and above inside K ( what they call 
the single coating) in Theorem 5, where, among other results, they prove that the tangential 
components of the electric and magnetic fields of their solutions have to vanish at dK^. They 
conclude that their solutions do not exists for generic currents inside and that cloaking 
holds for passive devices but that it fails for active devices with generic currents inside K. 
To deal with this issue they propose to add a perfect electric conducting lining to K -what 
makes it to appear as passive- or to introduce a different construction that they call the 
double coating. 

In this paper, as well as in ^4^ i5j , we proceed in a completely different way. Instead of 
postulating a priori a particular class of weak solutions in distributions sense across OK 
we first characterize all possible ways to define solutions that are compatible with energy 
conservation. They correspond to all self-adjoint extensions of the Maxwell generator. Note 
that each self-adjoint extension can be understood in terms of boundary conditions at dK±. 
We proved in [U [5] that all self-adjoint extensions are the direct sum of some self-adjoint 
extension inside K with a fixed self-adjoint extension outside K. This implies that the 
solutions inside and outside of K are completely decoupled from each other and that, in 
general, they are discontinuous at dK. Another consequence is that in the case where the 
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permittivity and the permeability are bounded below and above inside K weak solutions 
in distribution sense across dK, and more generally solutions with transmission conditions 
that link the inside and the outside of K, are not self-adjoint, i.e., they are not allowed by 
energy conservation. Note that in this the permittivity and the permeability are 

bounded below and above inside K, requiring that the tangential components of both, the 
electric and the magnetic field vanish at dK_ is not a self-adjoint boundary condition. We 
only are allowed to require that one of them vanishes. However, requiring that both vanish 
at dKj^ is a self-adjoint boundary condition because the permittivity and the permeability 
are degenerate at dK^. 

We also proved in [U [S] that cloaking of passive and active devices always holds for 
all possible ways to define solutions that satisfy energy conservation, i.e., with self-adjoint 
boundary conditions. 

There are many papers that discuss line transformed, or cylindrical, cloaks. See for 
example [IlElHEl[Zl[ini[IIl[I2l[Il[ISl[ISl[ISl[2ni. In [IB] boundary conditions are considered 
to enhance cloaking for a cylindrical cloak. Note that the results for line transformed -or 
cylindrical- cloaks are quite different from the ones for point transformed cloaks studied in 
this paper. 

As it is often the case in the papers on electromagnetic invisibility cloaks, we make the 
assumption that the media are not dispersive. This is a widely used idealization. As is 
well known, metamaterials are dispersive, and, furthermore, when the permittivity and the 
permeability have eigenvalues less than one, dispersion comes into play in order that the 
group velocity does not exceeds the speed of light. This idealization means that we have to 
take a narrow enough range of frequencies in order that we can analyze the cloaking effect 
without taking dispersion into account. In practice this means that cloaking will only be 
approximate. Note, moreover, that the results in this paper, as well as in [HIS], are proven 
for cloaks in exact transformation media. 

The paper is organized as follows. In Section 2 we prove that our cloaking boundary 
conditions at dK- are satisfied. In Section 3 we illustrate our method by considering one 
spherical cloak with an active device given by a radial electric current at the boundary of K. 
The case of a magnetic current at the boundary of K follows in the same way. We assume 
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that K is isotropic and spherically stratified. We verify in this particular case, by an explicit 
computation, that our cloaking boundary conditions are satisfied and that cloaking of active 
devices holds, even if the current is at the boundary of the cloaked object, as we have proven 
in Section 2 in the general case where there is no explicit solution. We end the paper with 
conclusions. 

For other results on invisibility cloaks see references [131 [ISl [Hj] as well as the references 
quoted there and in [3j . 

2 The Boundary Conditions 

Let us consider Maxwell equations in R^, in the time domain, 



VxE = -^B, VxH=-D. (2,1) 

V-B = 0,V-D = 0, (2.2) 

and in the frequency domain, assuming a periodic time dependence of E, H given by 
e™*, with u) the frequency, 

V X E = -iujB, V X H = iooB, ^ 0, (2.3) 
V-B = 0,V-D = 0, (2.4) 

where we have suppressed the factor e*'^* in both sides. In this paper we take the time factor 
e*"^* to use the convention of [3 [9]. Note that ([231) follows from fl23D . 

We briefiy recall some notations and definitions from [5]. 

Let us first consider the case where there is only one cloak located at x = 0. See Figure [TJ 
We designate the Cartesian coordinates of x by x'^, A = 1, 2, 3. To define the transformation 
media we introduce another copy of M^, denoted by Mq. The points in Mq are denoted by 
y with coordinates y'^,\ = 1,2,3. We designate, x := x/|x|,y := y/|y|. Consider the 
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following transformation from Mq \ {0} to M^, 



x = x(y) = /(y):=^?(|y|)y. (2.5) 

In spherical coordinates this transformation changes the radial coordinate but leaves the 
angular coordinates constant, i.e., |x| = 5f(|y|),x = y. Given < a < 6 we wish that this 
transformation sends the punctuated ball < |y| < 6 onto the concentric shell a < |x| < b, 
that it is the identity for |y| > & and that it is one-to-one. Then, we assume that g satisfies 
the following conditions. 

DEFINITION 2.1. For any positive numbers a, b with < a < b, we say the g is a cloaking 
function in [0,6] if g{p) is twice continuously differentiable on [0,b],g{0) = a,g{b) = b, and 
g'{p) :=A^(p)>0,pG[0,6]. 

We define. 



X = x(y) = /(y) := ^(|y|)y, forO < |y| < 6, 
X = x(y) := y, for |y| > 6. 

With these conditions fl2.6p is a bijection from Rq \ {0} onto R"' \ -Ba(O), where we denote 



S,.(xo) := {x G : I X - xol < r} . (2.7) 

Moreover, it blows up the point onto the sphere |x| = a. It sends the punctuated ball 
< |y| < 6 onto the concentric shell a < |x| <b and it is the identity for |y| > 6. It is twice 
continuously differentiable away from the sphere |y| = 6, where it can have discontinuities 
in the derivatives depending on the values of the derivatives of g at b. 

In [3] the quadratic case 



g{p)= l-^+p(p-6) 



p + a 



with p G R was discussed in connection with a cylindrical cloak in an approximate trans- 
formation medium. In pQ the first order case g{p) = ^^p + a was considered. First order 
transformations were previously used in pT| [22] in the context of Calderon's inverse conduc- 
tivity problem. 
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The closed ball i^' := {x G M'^ : |x| < a} is the region that we wish to conceal, and we 
call it the cloaked object. The spherical shell a < |x| < 6 is the cloaking layer. The union of 
the cloaked object and the cloaking layer is the spherical cloak. The domain |x| > 6 is the 
exterior of the spherical cloak. 

We now put a finite number of spherical cloaks in different points in space in such a 
way that they do not intersect. See Figure [21 Let us take as centers of the cloaks points 
Cj G M^, j = 1, 2, ■ ■ ■ where N is the number of cloaks and Cj ^ Ci,j 1,1 < j,l < N. We 
take < aj < bj, and cloaking functions (/j that satisfy the conditions of Definition 12.11 for 
Oj, bj,j = 1,2,3 ■■■ N, and we define the following transformation from Rq \ {ci, C2, ■ ■ ■ , c^v} 
to R^. 

X = x(y) = /(y) := Cj + gj{\y - Cj\) f^j,y E Bb^{Cj), j = 1,2,--- ,N, 

(2.8) 

X = x(y) = /(y) := y,y G \ uf^,B,^{c,), 
where B[,^{cj) are balls in Mq. 

The cloaked objects that we wish to conceal are given by, 



Kj := {x G : |x - Cjl < Oj} , j = 1, 2, ■ ■ ■ , N. (2.9) 

The spherical shells < |x — Cj| < hj,j = 1,2, ■ ■ ■ , N are the cloaking layers. The spherical 
cloaks are the balls Bh.{cj) in M^. We denote by K the union of all the cloaked objects, 

K := uf^.K,. (2.10) 

The domain 

Mnuf=i5.,(c,) (2.11) 

is the exterior of the all the spherical cloaks. We assume that the spherical cloaks are at a 
positive distance of each other, 

min distance (^Bb.{cj), Bb^{ci)) > 0, j 7^ /, j, / = 1, 2, ■ ■ ■ , A^. 

Denote 

Qo := Rl \ {ci, C2, ■ ■ ■ , c;v}, n:=R^\ K. 
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Then, f l2.8p is a bijection from onto Q, and for j = 1, 2, ■ ■ ■ , it blows up the point 
Cj onto the sphere |x — Cj| = aj. It sends the punctuated ball < |y — Cj| < bj onto 
the shell aj < |x — Cj| < bj and it is the identity for y G Mg \ interior (u^;^i?5^(cj)). It 
is twice continuously differentiable away from the spheres |y — Cj| = bj, where it can have 
discontinuities in the derivatives depending on the values of the derivatives of Qj at bj. 

The elements of the Jacobian matrix are denoted by Ay, 

Al, e [Qq \ Uf^^dBb^{Cj)), and that it can have discontinuities on UjLidBi,^{cj) depend- 
ing on the derivatives of gj at bj. We designate by A'^ the elements of the Jacobian of the 
inverse bijection, y = y(x) = /^-"^(x), 

- (2-13) 

A\ G (^7 \ U^]^9i?b^. (Cj)) , and it can have discontinuities on U^^^dBh_.{cj) depending on 
the derivatives of Qj at bj. 

It follows from (12.81) that the transformation matrix fl2.12p is given by. 



= i^^A^ + - '-^) (y - c,)^(y - c,)^' 



yeBh^icj),l<j<N, (2-14) 



(2.15) 



The determinant is equal to, 

A(y) = 9',{\y - c,|) y e i?,,(c,), i < j < iv, 

A(y) = l,yGM3\uf=ii?,^(c,). 
Note that A diverges at the boundary of K. 

We take here the material interpretation and we consider our transformation as a bijection 
between two different spaces, Qq and Q. However, our transformation can be considered, 
as well, as a change of coordinates in f2o- These two points of view are mathematically 
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equivalent. This means that under our transformation Maxwell equations in and in f2 
have the same invariance that they have under change of coordinates in three-space. See, 
for example, Let us denote by Eq, Hq, Bq, Dq, ^q*^, /ig'^, respectively, the electric and 

magnetic fields, the magnetic induction, the electric displacement, and the permittivity and 
permeability of Qq. e^'^jfiQ", are positive Hermitian matrices that are constant in Qq. The 
electric field is a covariant vector that transforms as, 

i?,(x) = ylf(y)Eo,A'(y). (2.16) 

The magnetic field H is a covariant pseudo- vector, but as we only consider space transfor- 
mations with positive determinant, it also transforms as in f l2.16p . The magnetic induction 
B and the electric displacement D are contravariant vector densities of weight one that 
transform as 

B\^) = {A{y))-'Ai,{y)B^'{y), (2.17) 

with the same transformation for D. The permittivity and permeability are contravariant 
tensor densities of weight one that transform as, 

.^^(x) = {A{y))-'AUy)A:,{y)e^'^'{y), (2.18) 

with the same transformation for Maxwell equations fl2.m2.4l) are the same in both 
spaces fl and Qq. Let us denote by sxu, fi'\u,£o\u, fJ'OXu, respectively, the inverses of the 
corresponding permittivity and permeability. They are covariant tensor densities of weight 
minus one that transform as. 



e,,(x) = A(y)A^'(yX(y)5oAv(y), /iA.(x) = A(y)A^'(yX(y) /.oAV'(y)- (2-19) 
We have that 

det e^'' = A-^ det det fi^" = A'^ det (2.20) 

det Exu = A det Eqxu, det fixu = A det fioxu- (2.21) 
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The matrices e^'^, are degenerate at OK and the matrices e\,j, ^xv are singular at dK. 

As the e^^ and yU'^^ are degenerate at the boundary of the cloaked object we have 
to make precise what do we mean by a solution to Maxwell equations in neighborhood of 
dK. In other words, we have to specify the cloaking boundary conditions that the solutions 
have to satisfy on the outside and on the inside of dK. We solved this problem in [4, ^ by 
requiring that the fundamental principle of energy conservation be satisfied. That is to say, 
we obtained the appropriate boundary conditions by requiring that the solutions conserve 
energy. Note that as our media are loss-less energy has to be conserved. Furthermore, any 
energy loss of the incoming waves could a priori be detected, and, in consequence, energy 
conservation is essential for cloaking purposes. We briefly review the results of [U [S]. 

We first consider the problem in VL. We write Maxwell equations in Schrodinger form. For 
this purpose we denote by e and /i, respectively, the matrices with entries Sx^ and ^xv Recall 
that (V X E)'^ = s^'^f'-^Ep, where s^'^'^ is the permutation contravariant pseudo-density of 
weight —1 (see section 6 of chapter II of [23], where a different notation is used). 

We define the following formal differential operator, 

-(h)-<4Vx"e)' (-) 

Here, as usual, we denote, eV x H := exuC^ x H)'', and /iV x E = fixuC^ x E)''. 
Equation (12. ip is equivalent to, 

4(h)--(h). P'^^) 

and equation (12.31) is equivalent to 




(2.24) 



Note that since the matrices e, fi are singular at dQ the operator has coefficients that are 
singular at dQ. This is the reason why we have to be careful when defining the solutions. 

It is necessary to define equation (I2.22p in an appropriate linear subspace of the Hilbert 
space of all finite energy fields in Q. We designate by Hqe the Hilbert space of all measurable, 
C^— valued functions defined on fl that are square integrable with the weight e^'^ and the 
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scalar product 



(E<".E<^')„,:= / 

Jn 



(2.25) 



Moreover, we denote by Tinn the Hilbert space of all measurable, C'^— valued functions 
defined on fl that are square integrable with the weight /i'^'^ and the scalar product, 



We first define in a nice set of functions where it makes sense, that we take as Cl{Q). 
In physical terms this means that we start with the minimal assumption that Maxwell's 
equations are satisfied in classical sense away from the boundary of Q. with domain 
D{aQ) := Cl{fl) is a symmetric operator in Hn, i-e. an C a^. To construct a unitary 
dynamics that preserves energy we have to analyze the self-adjoint extensions of a^, what in 
physical terms means that we have to make precise in what sense Maxwell's equations are 
solved up to dQ. In other words, to construct finite-energy solutions of (12.231) . with constant 
energy we have to demand that the initial finite energy fields , (E(0), H(O))"^ belong to the 
domain of one of the self-adjoint extensions of oq. The key issue is that has only one 
self-adjoint extension (i.e. it is essentially self-adjoint) that we denote by An. Moreover, An 
is unitarily equivalent to the free Maxwell propagator, Aq, in f^o- The unitary equivalence is 
generated by (I2.16P and by the same transformation for the magnetic field. This means that 
there is only one dynamics in Q that preserves energy, and that this dynamics is generated by 
An- As An and Aq are unitarily equivalent, the dynamics that they generate are physically 
equivalent, and this is the deep reason, from the point of view of fundamental physics, why 
there is perfect cloaking of passive and active devices. 

Solutions to (12.31 12.41) in general do not have finite energy because they do not have 
enough decay at infinity to be square integrable over all Q. Then, we only require that 
they are of locally finite energy in the sense that the electric and the magnetic fields are 
square integrable over every bounded subset of Q, respectively, with the weight e^"^, and 
/i^'^. Moreover, in order that the problem (12. 3[ 12.41) is well-posed -in the sense that it is 




(2.26) 



The Hilbert space of finite energy fields in Q is the direct sum 



'Hn ■= Ti-nE © T^qh- 



(2.27) 
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self-adjoint- the solutions with locally finite energy have to be locally in the domain of the 
only self-adjoint extension of a^^, that is to say, they have to be in the domain of Aq when 
multiplied by any continuously differentiable function with support in a bounded subset of 

n. 

On the basis of these considerations we proved in [H [5] that the solutions with locally 
finite energy in Q are solutions in distribution sense to (12.31 12.41) that satisfy, 

[ Exe^''Kdx^+ [ H X fi^" H^dx.^ < oo, (2.28) 
Jo Jo 

for every bounded set O G Q. Moreover, they have to satisfy the cloaking boundary condition, 

E X n = 0, H X n = 0, in dQ = dK+, (2.29) 

where dK+ is the outside of the boundary of the cloaked object and n is the normal 
vector to dK^. 

Note that as Aq is the only self-adjoint extension of a^, this is the only possible self- 
adjoint boundary condition on dK^. It is self-adjoint because the matrices e,fi are singular 
at dK+. Hence, cloaking as boundary value problem consists of finding a solution to (12. 3[ 
12. 4p in Q with locally-finite energy that satisfies the cloaking boundary condition given in 

□ 

Let us now consider the propagation of electromagnetic waves inside the cloaked object. 
We assume that in each Kj the permittivity and the permeability are given by Ej'^, i with 
inverses SjXui^jXv and where Sj^^j are the matrices with entries ejxu, fJ'jXu- Furthermore, we 
assume that < £:j"'(x), /ij"^(x) < C, x G Kj and that for any compact set Q contained in 
the interior of Kj there is a positive constant Cq such that det£:j''^(x) > Cg , det /ij''^ (x) > 
Cq, X G Q, j = 1, 2, ■ ■ ■ , A^. In other words, we only allow for possible singularities of ej, fij 
on the boundary of Kj. 

We designate by Hje the Hilbert space of all measurable, C^— valued functions defined 
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on Kj that are square integrable with the weight e^^ and the scalar product, 

(ES".Ef),,-X/S's"^<'-'- (2-30) 

Similarly, we denote by 'HjH the Hilbert space of all measurable, C'^— valued functions 
defined on Kj that are square integrable with the weight and the scalar product, 

K'.Hf)^^:= /^^ffj^f^rfx'. (2.31) 

The Hilbert space of finite energy fields in Kj is the direct sum 

Tij := 'HjE ® 'HjH, (2.32) 
and the Hilbert space of finite energy fields in the cloaked object K is the direct sum. 



Hk '■= ®j=iHj. 

The complete Hilbert space of finite energy fields including the cloaked object is. 



n:=nn®nK- (2.33) 

We now write (12. ip as a Schrodinger equation in each Kj as before. We define the 
following formal differential operator, 

(2.34) 



Equation (12. ip in Kj is equivalent to 

o o 

Let us denote the interior of Kj by Kj-= Kj \ dKj. Then, aj with domain Cl{Kj) is a 
symmetric operator in Tij. We denote, 

a := an © ax, where ax '■= ©^la^, (2.36) 
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with domain, 

D{a) { ( ) ®U ( H- ) ^ ®'^=' ^°^^^'^} • ^^-^^^ 

The operator a is symmetric in Ti. The possible unitary dynamics that preserve energy for 
the whole system, including the cloaked object, are given by the self-adjoint extensions 
of a. We proved in [U [5] that every self-adjoint extension, A, of a is the direct sum of Aq, 
and of some self-adjoint extension, A^-, of ax-, i.e., 

A = An® Ak. (2.38) 

This result implies that the cloaked object K and the exterior Vt are completely decoupled. 
That electromagnetic waves outside K can not go inside and vice versa that waves inside 
can not propagate outside, and that there is perfect cloaking of passive and active devices. 
Choosing a particular self-adjoint extension Ak amounts to fixing a boundary condition 
in the inside of the boundary of the cloked object, dK^. The self-adjoint extension -or 
boundary condition- that nature will take depends on the properties of the media inside the 
cloaked object. Note that this does not mean that we have to put any physical surface, a 
lining, on the surface of the cloaked object to enforce any particular boundary condition on 
the inside, since this plays no role in the cloaking outside. It is, however, of importance to 
determine what the self-adjoint extension in or the interior boundary condition, has to 
be for specific cloaked objects. See [HE] for a detailed discussion of these issues. 

The problem that we address on this paper is to determine what the boundary conditions 
in dK^ have to be in the particular case where the permittivity and the permeability in K 
are bounded and non-degenerate, i.e., when the matrices e^^ , are bounded above and 
below in Kj, 

< Ci < ef, /if < C2, X e i^j, J = 1, 2, ■ ■ ■ , N, (2.39) 

for some positive constants Ci,C2. This is clearly the most important case in the applica- 
tions. It corresponds to a standard object that is cloaked with a metamaterial. 

Let us consider the case of an active device with electric and magnetic currents in K. 
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The Maxwell equations at frequency uj are, 



V X H = iujT> + J, (2.40) 

V X E = -iu'B - Jm, (2.41) 

where J and Jm are, respectively, the electric and the magnetic currents, that we assume 
are different from zero only in K. 

As we mentioned above, we have already proven in [H[5] that energy conservation implies 
that the electromagnetic waves inside K can not propagate outside and that, vice versa, the 
waves outside can not go inside. The key issue here is that this is consistent with Maxwell 
equation (12.401) only if the normal component of the total current (i.e. the sum of the 
displacement current and the electric current) vanishes at dK^, i.e., if 

(icuD + J) ■ n\Q^^ = 0, (2.42) 

where as usual by dK^ we mean that we approach the boundary from the inside. In a similar 
way, the consistency with Maxwell equation (12.411) implies that, 

(icuB + J„) ■ n\gj^_ = 0. (2.43) 

Note that we do not need to ask that (12. 42^ I2.43P hold at dKj^^ because as we assume that 
J, Jm are identically zero outside K conditions (12. 42^ 12.431) in dK^ follow from equations 
(12.141 12.151 I2.17P and the same transformation equation for D, since the solution in Q is 
obtained applying the transformation formulae to a solution in Mg [HE]. Moreover, the 
boundary conditions (12.421 [2^43]) on dK_ and Maxwell equations (I2.40[ [2^411) imply that, 

(V X H) ■ n\gj^_ = 0, (V X E) ■ n\9j^_ = 0. (2.44) 

Hence, we have proven that the boundary conditions that we have to impose on the inside 
of the boundary of the cloaked object are (I2.44p . namely that the normal components of the 
curl of the electric and the magnetic fields have to vanish. The Maxwell propagator ax with 
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the boundary condition fl2.44l) has aheady been studied in the mathematical hterature, in 
particular in relation with Beltrami fields. As it was to be expected from our analysis, 
with the boundary condition (12.441) is a self-adjoint operator. In other words, the boundary 
condition fl2.44l) defines the self-adjoint realization, A^-, of the Maxwell propagator in K 
that is imposed by energy conservation. For the proof of self-adjointness, as well as other 
issues, including the formulation of fl2.44p in weak sense see [HI [251 EZj, in particular, 
see page 158, Theorem 2.1, Corollary 2.1.1, page 164, and Theorem 2.3 of [26j. Note that in 
our case the Neumann fields are zero. Remark that imposing that the normal components 
of D and B are zero at dK^ is not a self-adjoint boundary condition, i.e., it does not define 
a self-adjoint extension of Maxwell generator in K. 

We have now a complete formulation of cloaking as a boundary value problem. It consists 
of finding a solution of Maxwell equations (12.401 12.4ip in distribution sense in \ dK, with 
locally finite energy, i.e., they satisfy 

! Exe^''Kdx.^+ [ H X fi^"" H^dx.^ < oo, (2.45) 
Jo Jo 

where O is any bounded subset of R^. Moreover, they have to satisfy the cloaking boundary 
conditions 

E X n = 0, H X n = 0, at dQ = dK+, (2.46) 

and 

(V X E) • n = 0, (V X H) • n = 0, atdK_. (2.47) 

We have derived the boundary conditions fl2.46[ I2.47P by requiring that the solutions 
to the fixed frequency Maxwell equations (12.401 I2.4ip are (locally) in the domain of the 
appropriate self-adjoint extension (12.380 of Maxwell generator fl2.36p . Note that when we 
define the self-adjoint operator Ax we have to require that all functions on its domain 
satisfy the boundary conditions (12.470 . not just the solutions to the fixed frequency Maxwell 
equations. In fact, the choice of the self-adjoint Maxwell generator Ax has implications that 
go well beyond the formulation of cloaking as a boundary value problem. For example, it 
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determines the time evolution of finite-energy wave packets in the time domain. See [U [3] 
for this issue. 

3 The Case of a Radial Source at the Boundary of K 

In this section we illustrate our method by considering an active device given by a radial 
electric current at the boundary of K. The case of a magnetic current at the boundary of K 
follows in the same way. We assume that K is isotropic and spherically stratified, i.e., that 
the permittivity and the permeability depend only on |x|. The case where K is isotropic 
and homogeneous, and with an electric dipole contained in the interior of K was already 
considered in [14|. We verify in this particular case, by an explicit computation, that our 
cloaking boundary conditions are satisfied and that cloaking of active devices holds, even if 
the current is at the boundary of the cloaked object, as we have proven in Section 2 in the 
general case where there is no explicit solution. 

We assume that we have only one spherical cloak, located at the origin, i.e., = 
l,Ci = 0, that K is isotropic with permittivity and permeability, that are bounded, 

that they have a positive lower bound, and that they depend only on r := |x|. For simplicity 
we take a first order transformation with g{p) = + a,0 < a < b. In the cloaking layer 
the permittivity and the permeability tensors are given by [H, [16] . 

e^" = ErX.^x.'' + Ete^r + = /u.x^x'^ + ptO^t + /it<^^^^ a < |x| < 6, (3.1) 

where, 9, (p are unit tangent vectors, respectively, to the coordinate lines, r, constant and 
r, Q constant, in spherical coordinates, r, 6', Moreover, 

Etje^ = fit/fJ'O = b/(b- a),er/et = Pr/fJ't = (r - af /r^,a < r < b. (3.2) 

We assume that for r > b the medium is homogeneous and isotropic with permittivity and 
permeability, Eq, Pq. 

The expression of the transverse electric, TE, and transverse magnetic, TM, fields in 
terms of potentials given in Section 8.6 of [S] remain true in our case (remark the Et and /^^ 
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are constant in the cloaking layer). See also [13], [15]. TE and TM fields decouple, and since 

we have a radial electric current we only consider TM modes. Assuming that J = Jr{r)'k 

and that = 0, the TM fields are given by the potential as follows [H], 

—i ^ d 

Et = gradg — /, 

uJEtr ^ dr 

H* = ^ ( grades / x x) , 

= —- -IS.Qu}I — : Jr- (3-3) 

We expand the potential in spherical harmonics, 

mn 

and we assume that the radial current has the following expansion, 

Jr = 5Zj,nn5(r-a)Kr(^,¥'), 
mn 

for some constants Jmn- For example, for an electric dipole located at (0, 0, a) we have that 

J _ iu)Pe 1 2n+l X rcl 
■Jmn— 27ra2'\/ 47r ""^.O l^" 

We now set the inner boundary at a + 5, for small 5 > 0, i.e., we assume that the 
permittivity and permeability are equal to for r < a + 5, that they are given by (13. 

13. 2p ioi a + 5 < r < h and by Eq, for r > b. We compute the solution, and then we take 
the limit as 6 tends to zero ([19]). 



For < r < a+6, the potential Ia,mn created by the source satisfies the following equation 
A 

dr 

Ia,mn Jmn^iX '^)- (3-4) 



[8], where we denote -r; by ' 



1 , 2 + 1) 

U^Ellil ^ 



T" - —e' r + 

a,mn 1 a,mn ' 

^1 



Let Vn and Wn be to independent solutions of the homogeneous equation with f„ regular at 
zero. For example, if £i(r),/ii(r) are piecewise constant, f„ can be taken as a Ricatti-Bessel 
function of the first kind and Wn as a Ricatti-Bessel function of the second or third kind [9] 
in each layer where Ei,fii are constant. Hence, the solution to (13. 4p is given by. 



w„(a)f„(r),0 < r < a, 
Vn{a)Wn{r), a < r < a + 6. 



v'^[a)w„[a) - Vn[a)w'^[a] 



19 



We suppose that Vn{a) and m„(a) are different from zero. We assume that there is also 
a reflected wave. Then, the total potential is given by, 



(r) + RmnVnir), < r < a + 6, 



(3.6) 



where the Rmn are the reflection coefficients. In the cloaking layer the potential satisfies the 
equation (remark that et, fj^t are constant). 



2 St nin + 1) 



0,a + 6 <r <b. 



Br 

The solution is given by Ricatti-Bessel functions of the first and second kind, 



(3.7) 



Imnir) = Cmni'niktir - tt)) + dmnXnihir - o)), h := uj^/s^t, a + 6 < r < b. (3.8) 

Outside of the cloaking layer the potential satisfies equation (13.71) with et = Sr = Eq and 
fit = /io- The solution is an outgoing wave, 

Imn{r) = Tmn (n{kor) , /Cq := Uy/EofiQ, b < r, (3.9) 

where C is a Ricatti-Bessel function of the third kind. 

Requiring that the tangential components of the electric and the magnetic fields are 
continuous at a + 6 and b we obtain the following equations. 



1 , kt 

' I'^„^^{a + 6) + RmnV'nia + 6)) = — (Cmn^niktS) + dmnX'niktS)) , (3.10) 



^ /" I X\ \ a,mn\ ' ' -^"inii^ny ' 'JJ ^ 

ei[a + d) ' €t 

h,mn{a + 5)+ RmnVn{a + 5) = C.mnAi{ktS) + dmnXniktS) , (3.11) 

k k 

— {Cmn^nikt{b - a)) + dmnX'niktib - a)) = —TmnCnikob) , (3.12) 

Cmni^n{kt{b - a)) + dmnXn{kt{b - a)) = TmnCn{kob) . (3.13) 

Let us denote, 

a := I 



a:={^y'\ (3.14) 



A := Cikob)Xn{kt{b - a)) - aCn{kob)Xn{kt{b - a)), (3.15) 
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132 ■= C'ni.koh)'ilJn{kt{h - a)) - aCnikob)ip'^{kt{b - a)), := (3.16) 

/?4 := <(a + 6mMhS) + Xnih6)) - + _^ 5)(/33^;(M) + XnihS)), (3.17) 

Solving (13.10113. 15]) in the generic case where 7^ we obtain that, 



Rmn = r— TTT iiPsi^nihS) + Xn{kt5))dmn - Ia,mn{a + 6)) , (3.18) 



Cmn (3^dmni (3.19) 



(a)-, (3.20) 

Tmn = -—^r-r:{f33'ipn{kt{b - a)) + Xnihib - a))d^n- (3.21) 

Using the expansions of the spherical Bessel functions for small argument given in equa- 
tions (10.1.2, 10.1.3) of [2S] we prove that, 

dmn = -Jmn , , ^ ""f""]'' TTTT (^*^)"^' + 0^5)-^'). (3.22) 

ktei{a + o)Vn\a + d)n[2n — 1)!! 

Using (13.221) we obtain the small 5 expansions of Rmn,Cmn,Tmn, respectively, from f l3.18[ 
1339113:211) . 



Rmn= / , (-/a,n^n(a + S) + 0{kt6), (3.23) 

Vn[a + d) 

/33 ( -Jmn , , "ff' nn + OiiktSr'-')) , (3.24) 

ktei{a + o)Vn{a + d)n[2n — iy.\ J 



T 



JmnVn{o)et{l3?,i'n{kt{b - o)) + Xn{kt{b - O)) 



i;:ia:^;cMa:s):^:- m " (m)- (3.25) 
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Moreover, by ^ EUl 

/™„(a + S) = Jrnn , , ^"^'^^^^ , M + (3.26) 
■u„(a + d)nktei{a + d) 

It follows that when 5 = 0, Cmn = c^mn = = 0, and the electric and magnetic fields 
outside K are zero, as predicted by our theoretical results. Also, for 5 = 0, Imn{,o) = 0, and 
then, by equation (13.31) the radial component of the sum of the displacement current and 
the electric current vanishes at r = a. By Maxwell equation (12.401) the normal component 
of the curl of the magnetic field vanishes ar r = a. The corresponding statement for the 
normal component of the curl of the electric field is trivial in this case by Maxwell equation 
(12.411) and as the magnetic field is transversal and the magnetic current is zero. Hence, 
the boundary conditions (I2.44p are satisfied and cloaking holds even if the current is at the 
boundary of K, as predicted by our theoretical results. 



4 Conclusions 



The results of this paper and of [U [5] give a complete rigorous mathematical analysis of 
point transformed electromagnetic invisibility cloaks. We solved the mathematical chal- 
lenges posed by the fact that the permittivity and the permeability are degenerate at the 
boundary of the cloaked object K. In particular, we characterized all possible ways to define 
solutions of Maxwell equations that are compatible with energy conservation. This result 
was obtained by characterizing all possible boundary conditions at dK± that are allowed 
by energy conservation. They correspond to all self-adjoint extensions of the Maxwell gen- 
erator. As it turned out, all self-adjoint extensions are the direct sum of some self-adjoint 
extension inside K with a fixed self-adjoint extension outside K. This implies that the 
solutions inside and outside of K are completely decoupled from each other and that, in 
general, they are discontinuous at dK. We also proved that cloaking of passive and active 
devices always holds for all possible ways to define solutions that satisfy energy conservation, 
i.e., with self-adjoint boundary conditions. The boundary condition at dK^ is always that 
the tangential components of both the electric and the magnetic field vanish. At dK_ the 
boundary condition can by any self-adjoint boundary condition for the Maxwell generator 
in K. The particular self-adjoint boundary condition that nature will take depends on the 
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specific properties of tlie media inside K. In tliis paper we solved tfie problem of determin- 
ing the appropriate boundary conditions at dK_ in the case where the permittivity and the 
permeabihty inside K are bounded and have a positive lower bound. This corresponds to the 
situation where we have a standard object K -that could be anisotropic and inhomogeneous, 
but whose permittivity and permeability are neither singular nor degenerate- that is coated 
by a transformation medium that is degenerate at dK^ . We also allow for passive and active 
devices in K. This is perhaps the more important case in the applications. 

In this way we have obtained a complete formulation of cloaking with passive and active 
devices as a boundary value problem. It consists of finding a solution of Maxwell equations 
fl2.40l 12.411) at frequency u, in distribution sense in \ dK, with locally finite energy, i.e., 
they satisfy 

[ Exe^^E^dx^ + [ Hxfi^'^H^dx^ < oo, 
Jo Jo 

where O is any bounded subset of M^. Moreover, they have to satisfy the following cloaking 
boundary conditions 

E X n = 0, H X n = 0, at dK+, (4.1) 

and 

(V X E) ■ n = 0, (V X H) ■ n = 0, &tdK^. (4.2) 

In the case of one spherical cloaked object, that is isotropic and spherically stratified, and 
that has an active device given by a radial electric current at the boundary, we verified by an 
explicit computation that our cloaking boundary conditions are satisfied and that cloaking 
of active devices holds, even if the current is at the boundary of the cloaked object. 

A novel aspect of our work is that we proved our results for transformations media that 
are obtained from general anisotropic media, i.e., that it is not necessary to transform from 
isotropic media. This means that it is possible to cloak objects that are contained inside 
general anisotropic materials, general crystals for example. For this purpose, we just have 
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to take as the permittivity and the permeabihty of the general anisotropic medium before 
transformation the ones of the general anisotropic material, or the general crystal, that 
contains the object that we wish to cloak. The fact that it is possible to cloak objects 
inside general anisotropic media opens the way to other interesting potential applications, 
for example to guide electromagnetic waves under quite general circumstances. 

The results above -that are proven for cloaks in exact transformation media (ideal cloaks)- 
set the stage for the rigorous study of the cloaks in the approximate transformation media 
that one has to consider in practical situations, what is one of the main open questions on 
this area. This issue can be understood as the problem of the stability of cloaking under 
the perturbation on the permittivity and the permeability given by the difference in the 
permittivity and the permeability between the exact and the approximate transformation 
media. Our formulation of cloaking as a self-adjoint problem shows that the important issue 
of stability of cloaking can be formulated as a problem in perturbation theory of the self- 
adjoint Maxwell generator. Perturbation theory of self-adjoint operators is a main stream 
topic in modern mathematical physics and there is a large body of results. See for example 
[29] . Our analysis opens the way to a rigorous study of the stability of cloaking along these 
lines. 

For exact transformation media the boundary conditions (14.11 14. 2p are satisfied because 
they follow from energy conservation and there is no need to add any lining to impose them. 
However, our results suggest a method to enhance cloaking in the approximate transfor- 
mation media that are used in practice. Namely, to coat OK by a material that imposes 
the boundary conditions (14. ![ 14. 2p . As these boundary conditions have to be satisfied for 
exact transformation media, adding a lining that enforces them in the case of approximate 
transformation media will improve the performance of approximate cloaks. 
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